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ABSTRACT 

We report the results of the narrow-band photometry and imaging monitoring of comet C/2013 R1 (Lovejoy) with the robotic telescope 
TRAPPIST (La Silla observatory). We gathered around 400 images over 8 months pre- and post-perihelion between September 12, 
2013 and July 6, 2014. We followed the evolution of the OH, NH, CN. C3, and C2 production rates computed with the Haser model, as 
well as the evolution of the dust production. All five gas species display an asymmetry about perihelion, since the rate of brightening 
is steeper than the rate of fading. The study of the coma morphology reveals gas and dust jets that indicate one or several active 
zone(s) on the nucleus. The dust, C2, and C3 morphologies present some similarities, while the CN morphology is different. OH and 
NH are enhanced in the tail direction. The study of the evolution of the comet activity shows that the OH, NH. and C 2 production rate 
evolution with the heliocentric distance is correlated to the dust evolution. The CN and, to a lesser extent, the C3 do not display such 
a correlation with the dust. This evidence and the comparison with parent species production rates indicate that C2 and C3, on one 
hand, and OH and NH, on the other, could be -at least partially- released from organic-rich grains and icy grains. On the contrary, all 
evidences point to HCN being the main parent of CN in this comet. 
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1. Introduction 

Comets are remnants of the early solar system. They have un¬ 
dergone little alteration since their formation, and their nuclei 
are thought to preserve physical and chemical information about 
the early solar nebula. When a comet approaches the Sun, the 
ices contained in the nucleus sublimate, and complex molecules 
are released into the coma. These molecules, the parent species, 
can be observed in the infrared and radio domains. In optical 
studies we observe daughter or even granddaughter species (CN, 
C 3 , C 2 , etc.) that are the product of the parent species’ photodis¬ 
sociation. Since the relationship between parent, daughter, and 
granddaughter volatiles is still not understood well, it remains 
very difficult to link the radicals observed in the optical to nu¬ 
cleus ice abundances. Recent studies comparing the abundances 
of parent volatiles derived from infrared observations and abun¬ 
dances of daughter species derived from optical spectroscopy or 
photometry reveal a certain diversity among comets (Mumma & 
Charnley 2011). 

The case of CN is a good illustration of the problems en¬ 
countered while trying to link parent and daughter species. CN 
was first thought to be the product of HCN photodissociation. 
However, studies have shown that for some comets, the HCN 
abundance was insufficient to account for the CN production 
rates (Fray et al. 2005, Dello Russo et al. 2009). The same dis¬ 
crepancies were also noted between C 2 H 2 and C 2 production 
rates. Moreover, the spatial distribution of C 2 and CN in the 
coma is sometimes not consistent with these two species be¬ 
ing the product of C 2 H 2 and HCN photodissociation. It has then 
been suggested that these radicals could have other parents or 
could be released by organic-rich grains in the coma (A’Hearn 


et al. 1986, Combi & Fink 1997). However, additional parents 
for these two species are not needed for all comets. Since the 
sample of comets observed simultaneously in the infrared and 
the optical is still scarce, it remains difficult to have a compre¬ 
hensive picture of the relationship between parent and daughter 
species. 

We present here the monitoring of comet C/2013 R1 
(Lovejoy) performed with TRAPPIST (Jehin et al. 2011). We 
followed the evolution of the production rates of five daughter 
species and their coma morphology. We explore the link between 
the comet activity evolution and morphology and the possible 
origin of the observed radicals. Comet C/2013 R1 (Lovejoy) is 
a dynamically old long-period comet discovered on September 
7, 2013 at 1.97 au from the Sun (Guido et al. 2013). It was first 
detected by the Australian amateur astronomer Terry Lovejoy at 
magnitude 14.4 and brightened up to magnitude 4.5. The comet 
reached perihelion on December 22, 2013 at 0.8 au from the 
Sun. It follows an inclined and highly eccentric orbit (i = 64° 
and e = 0.9924). After the disintegration of comet ISON in late 
November, comet Lovejoy became an ideal backup target, and 
it has been extensively observed with a variety of instruments at 
different wavelengths around its perihelion passage. Thanks to 
the fast response possible with TRAPPIST, we were able to al¬ 
ready observe the comet a few days after its discovery. Since the 
comet was also observable post-perihelion, it was an extremely 
interesting target for the study of its activity and chemical com¬ 
position, and we started an intensive monitoring with narrow- 
band filters. 
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2. Observations and data reduction 

We observed comet Lovejoy with the TRAPPIST 60-cm robotic 
telescope at La Silla observatory (Jehin et al. 2011). TRAPPIST 
is equipped with a 2Kx2K thermoelectric ally cooled FLI Proline 
CCD camera with a field of view of 22’X 22’. We binned the 
pixels 2 by 2 and obtained a resulting plate scale of 1.302”/pixel. 
We observed the comet with HB narrow band filters (Farnham 
et al. 2000) isolating the emission bands of OH, NH, CN, C 3 , 
and C 2 , as well as emission free continuum at four wavelengths 
(see Farnham et al. 2000 for filters characteristics). Images have 
also been taken with broad band B, V, Rc, and Ic Johnson-Cousin 
filters. 

We collected about 400 images over eight months. We 
started the observations on September 12, 2013 when the comet 
was at 1.91 au from the Sun. We observed once or twice a week 
until November 16, when it was too close to the Sun. Just after 
the perihelion, the comet was too far north to be observed from 
La Silla, and we recovered it on February 13, 2014. We kept ob¬ 
serving with narrow band filters as long as the comet was bright 
enough, until July 6 , 2014 (r — 3.10 au, r being the heliocen¬ 
tric distance). After April 25 (r = 2.20 au) the signal was too 
low in the OH and NH filters to derive reliable production rates. 
Exposure times ranged from 15 to 120 seconds for broad band 
filters and from 60 to 900 seconds for narrow band filters. Most 
nights were photometric, and we discarded cloudy nights. 

Table 1 . Parent ( L P ) and daughter (L D ) scalelengths used in the 
Haser model. 


Species 

Lp 

(km) 

Ld 

(km)) 

OH 

2.40 10 4 

1.60 HP 

NH 

5.00 10 4 

1.50 10 5 

CN 

1.30 10 4 

2.10 10 5 

c 3 

2.80 10 3 

2.70 10 4 

C 2 

2.20 10 4 

6.60 10 4 


Notes. These values are given for 1 au and scale as r 2 (A’Hearn et al. 
1995). 


Calibration followed standard procedures using frequently 
updated master bias and flat and dark frames. The removal of the 
sky contribution and the photometric calibration of the images 
followed the procedure described in Opitom et al. (2015). After 
calibration, we derived a median radial profile for each image. 
We then removed the dust contamination in the gas radial pro¬ 
files by subtracting a properly scaled dust profile. From the OH, 
NH, CN, C 3 , and C 2 radial profiles, we derived production rates 
by adjusting a radial coma brightness model. The model adjust¬ 
ment is performed at a nucleocentric distance around 10 000 km 
to avoid PSF and seeing effects around the optocenter. At greater 
nucleocentric distances, the signal usually becomes fainter, es¬ 
pecially in the OH filter for which the signal-to-noise ratio 
is lower. Fluorescence efficiencies (also called g factors) from 
David Schleicher’s website 1 were used to convert fluxes into 
column densities. C 2 g factors only consider the C 2 in triplet 
state (A’Hearn 1982). C 3 g-factors are also from A’Hearn (1982). 
CN and NH fluorescence efficiencies vary with both the helio¬ 
centric distance and velocity and are respectively taken from 
Schleicher (2010) and from Meier et al. (1998). For OH g fac- 

1 Link: http://asteroid.lowell.edu/comet/gfactor.html 


tors, which also vary with the heliocentric velocity, we con¬ 
sidered unquenched values of the ground state (Schleicher & 
A’Hearn 1988). All fluorescence efficiencies are scaled as r 2 . We 
used the Haser model to compute gas production rates (Haser 
1957). This model is fairly simple and considers the photodis¬ 
sociation of mother molecules into daughter and granddaugh¬ 
ter molecules in a spherically symmetric coma. The model has 
only two free parameters: the mother and daughter scalelengths. 
We used the scalelengths from A’Hearn et al. (1995). These are 
shown in Table 1 for a comet at 1 au from the Sun and are scaled 
as r 2 . We used a constant velocity of 1 krn/s for both parent and 
daughter species, as in A’Hearn et al. (1995). 

From the profiles in the dust filters, we derived the A(9)fp 
parameter at three wavelengths. This quantity, first introduced 
by A’Hearn et al. (1984), is representative of the dust activity 
and independent of the aperture size if the dust coma is in a 
steady state (following a 1 /p radial profile, p being the nucleo¬ 
centric distance). Since the phase angle was varying significantly 
during the observations (from almost 70° around mid-November 
2013 to only a few degrees in early July 2014), we corrected the 
A(9)fp from the phase angle effect to obtain A(0)fp. Several 
phase functions have been successively proposed to correct 
the A(9)fp for the phase angle effect (Divine 1981, Hanner & 
Newburn 1989, Schleicher et al. 1998 or Marcus 2007). We used 
the phase function described by David Schleicher 2 , which is rep¬ 
resented in Fig. 1. It is a composite of two different phase func¬ 
tions from Schleicher et al. (1998) and Marcus (2007). We nor¬ 
malized the A(9)fp values at a 0° angle. 



Phase Angle [°] 

Fig. 1. Phase function used to correct the A(9)fp from the phase 
angle effect. 


3. Data analysis and results 

In this section, we present the evolution of comet Lovejoy activ¬ 
ity and composition with heliocentric distance. We compare the 
activity before and after the perihelion passage and also discuss 
the morphology of the dust and gas coma. 

2 Link: http://asteroid.lowell.edu/comet/dustphase.html 
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Table 2. OH, NH, CN, C 3 , and C 2 production rates and Afp for comet C/2013 R1 (Lovejoy) 


UT date 

r 

(au) 

A 

(au) 

Q(OH) 

Production rates (10 23 
Q(NH) Q(CN) 

mol/s) 

Q(C 3 ) 

Q(C 2 ) 

A(0)fp 
(10 2 cm) 

2013 Sep 12.38 

1.91 

1.81 

870±135 



0.67±0.04 

2.35±0.13 


2013 Sep 16.39 

1.85 

1.71 



3.14±0.12 


2.84±0.20 

3.79±0.28 

2013 Sep 21.33 

1.79 

1.59 

823±144 

4.89±1.63 

3.17±0.14 

0.74±0.06 

3.10±0.21 

4.37±0.37 

2013 Sep 24.39 

1.75 

1.52 



3.67±0.12 




2013 Sep 25.35 

1.73 

1.50 



3.70±0.14 


4.32±0.15 

4.04±0.39 

2013 Oct 04.35 

1.61 

1.28 


9.73±1.13 

3.93±0.12 




2013 Oct 05.33 

1.60 

1.25 

1120±180 


4.28±0.13 

1.03±0.04 

5.29±0.15 

6.30±0.18 

2013 Oct 08.33 

1.56 

1.18 

1520±180 


3.74±0.12 


4.49±0.16 

6.85±0.24 

2013 Oct 10.31 

1.53 

1.13 

1470±210 


5.02±0.16 


5.26±0.42 

7.52±0.24 

2013 Oct 20.29 

1.40 

0.89 

2160±330 


5.67±0.23 


7.45±0.34 

10.00±0.50 

2013 Oct 21.33 

1.38 

0.87 

2380±250 

15.80±1.70 

6.38±0.19 

1.19±0.04 

7.63±0.21 


2013 Oct 24.32 

1.34 

0.80 

2720±290 

19.10± 1.80 

6.96±0.21 

1.85±0.08 

9.10±0.34 

13.10±0.30 

2013 Oct 26.30 

1.32 

0.75 

2860±330 

20 .00±2.00 

7.20±0.31 

1.93±0.07 

9.62±0.39 

14.20±0.40 

2013 Nov 02.35 

1.23 

0.61 



7.97±0.27 

1.92±0.07 

12.20±0.30 


2013 Nov 03.37 

1.21 

0.59 

3010±350 

23.70±1.70 





2013 Nov 13.34 

1.09 

0.43 

3240±790 

31.10±4.60 

10.10±0.70 


16.10±0.50 

22.60±0.80 

2013 Nov 16.34 

1.06 

0.41 


43.80±8.60 

12.90±1.20 

3.40±0.31 

20.60±0.09 

28.50±2.10 

2014 Feb 13.39 

1.26 

1.56 



8.76±0.44 


13.70±0.50 

42.00±1.70 

2014 Feb 15.40 

1.28 

1.56 

1630±320 

20.80±3.50 



11.30±0.40 

34.90±1.40 

2014 Feb 16.40 

1.30 

1.57 



6.12±0.36 




2014 Feb 21.39 

1.36 

1.58 

1380±230 

17.50±2.60 

5.97±0.26 

2 .10±0.11 

8.63±0.31 

28.30±0.90 

2014 Feb 22.39 

1.37 

1.58 





9.59±0.27 

29.10±1.00 

2014 Feb 28.39 

1.45 

1.59 






21.70±0.90 

2014 Mar 01.39 

1.47 

1.59 






24.00±0.90 

2014 Mar 06.39 

1.54 

1.59 

1320±160 

14.30±1.70 

5.10±0.17 

1.69±0.06 

7.10±0.17 

21.80±0.60 

2014 Mar 13.39 

1.63 

1.58 

1160± 130 

11.90±1.50 

4.45±0.16 

1.46±0.06 

5.87±0.17 

17.80±0.60 

2014 Mar 17.39 

1.68 

1.57 






15.20±1.20 

2014 Mar 18.38 

1.70 

1.57 

939±115 

9.10±1.76 

3.99±0.17 

1.33±0.08 

5.97±0.26 

12.50±0.70 

2014 Mar 20.34 

1.72 

1.56 



4.01±0.14 




2014 Mar 25.39 

1.79 

1.55 



3.80±0.12 


4.66±0.21 

13.20±0.50 

2014 Mar 26.38 

1.80 

1.55 

859±91 

8.88±1.16 

3.95±0.14 

1.25±0.05 

4.43±0.29 

13.10±0.40 

2014 Apr 10.38 

2.00 

1.50 

627±68 


3.40±0.14 

0.98±0.04 

3.26±0.19 

9.35±0.35 

2014 Apr 14.42 

2.06 

1.49 



2.84±0.15 


2.72±0.23 

8.92±0.48 

2014 Apr 15.25 

2.07 

1.48 

561±56 



1.02±0.03 



2014 Apr 16.42 

2.08 

1.48 


2.93±2.45 





2014 Apr 18.38 

2.11 

1.48 

512±60 


2.72±0.15 


2.91±0.24 

7.56±0.45 

2014 Apr 19.34 

2.12 

1.47 


3.79±1.69 


0.72±0.07 



2014 Apr 22.30 

2.16 

1.47 

476±56 


2.93±0.13 


3.07±0.30 

7.21±0.40 

2014 Apr 25.38 

2.20 

1.47 

413±65 


2.50±0.13 


2.03±0.29 

6.70±0.33 

2014 May 14.16 

2.44 

1.50 



2.22±0.23 


1.72±0.40 

4.54±0.76 

2014 May 24.34 

2.57 

1.57 



1.80±0.21 



3.96±0.39 

2014 May 27.32 

2.61 

1.60 



1.72±0.18 

0.68±0.05 

1.97±0.45 

3.66±0.39 

2014 May 31.18 

2.66 

1.64 



1.91±0.14 

0.71±0.07 

1.98±0.30 


2014 Jun 07.21 

2.74 

1.73 






3.61±0.26 

2014 Jun 30.21 

3.02 

2.15 






3.01±0.41 

2014 Jul 04.03 

3.07 

2.23 






3.05±0.35 

2014 Jul 06.99 

3.10 

2.30 






2.62±0.59 


Notes, r and A are respectively the heliocentric and geocentric distances. The date given in the first column is the mid-time of the observations. 


3.1. Production rates evolution 

This section focuses on the analysis of the comet activity and 
composition evolution with the heliocentric distance. For each 
image in the narrow band filters, we computed the A(6)fp value 
or the production rate, depending on the filter. The results are 
given in Table 2, where we give the normalized A(0)fp value 
and the production rates for each image. Gas species are also 
represented all together in Fig. 2. The errors take both the errors 
in the flux calibration process and the sky contribution removal 
into account. Figure 3 shows the evolution of the OH, NH, CN, 
C 3 , and C 2 production rates and the A(0)fp as a function of the 
heliocentric distance along with the fitted slope for the evolution 
of pre- and post- perihelion activity. From both the gas and the 


dust production evolution, we see that the comet activity evolves 
quite regularly. We could not detect any major outburst or sudden 
change of the comet activity with the heliocentric distance. 

We immediately notice that the five gas species and the dust 
display an asymmetry about perihelion, since the rate of bright¬ 
ening is steeper than the rate of fading. The asymmetry does not 
have the same amplitude for the gas and the dust. Indeed, the 
A(0)fp values are higher after perihelion than before at the same 
heliocentric distance by at least a factor two, while the gas pro¬ 
duction rates do not display such a strong difference between 
pre- and post-perihelion values. This of course implies that the 
dust-to-gas ratio is higher after perihelion. 
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Fig. 2. Logarithm of the OH, NH, CN, C 3 , and C 2 production 
rates as a function of the heliocentric distance r (au). The /--axis 
is interrupted at perihelion ( 0.8 au). 


Table 3. Pre- and post-perihelion fitted power-law slopes for 
OH, NH, CN, C 3 , and C 2 production rates and A(0)fp helio¬ 
centric dependences. 


Species 

r-dependence 

pre-perihelion post-perihelion 

OH 

-3.23±0.33 

-2.62+0.19 

NH 

-3.55±0.28 

-2.75+0.31 

CN 

-2.60±0.17 

-1.84+0.09 

c 3 

-2.78±0.40 

-1.75+0.12 

c 2 

-3.56+0.16 

-2.76+0.12 

A(0)fp 

-3.92+0.19 

-3.01+0.07 


We adopted a power-law fit to represent the heliocentric de¬ 
pendence of the production rates and the A(0)fp both pre- and 
post-perihelion. The values of the slopes are given in Table 3. 
The rate of brightening is steeper than the rate of fading, which 
confirms the asymmetry about perihelion observed in Figure 3. 
The slopes have similar values for OH, NH, and C 2 and are con¬ 
sistent within the error bars with the slope observed for the dust 
heliocentric dependence. However, the slopes for CN and C 3 are 
different from the other gas species and the dust. This difference 
is visible on both sides of perihelion. 

We have a few observations of comet Lovejoy with the RC 
and GC continuum filters in addition to the ones regularly per¬ 
formed with the BC filter. The comparison of the A(6)fp values 
in these three filters allows the color of the dust and its evolution 
to be studied. The dust color is usually computed as the normal¬ 
ized gradient of A(6)fp from two continuum filters: 


color[/li, A 2 ] 


A fp\ ~ A fPi 2000 
Ti - A 2 Afpi+Afp 2 ' 


( 1 ) 


The normalized reflectivity gradient is expressed as the percent¬ 
age of reddening by 1000 A. We measured values between 6 % 


and 14 %/1000 A. This is in the range of what is usually ob¬ 
served for comets (see for example Lin et al. 2013; Lin et al. 
2012; Lara et al. 2004). The dust color stayed constant during 
our observations. 

3.2. Composition 

We monitored several production rate ratios to study the evolu¬ 
tion of the coma composition. In Fig. 4 we show the evolution 
of the C 2 -to-CN ratio. This ratio decreases with the heliocentric 
distance by almost a factor two between 1.1 and 2.5 au. At large 
heliocentric distances, the ratio does not seem to vary anymore, 
but the error bars are larger. Pre- and post-perihelion ratios match 
within the error bars, and both display the same trend. The C 2 
production rate is rising faster than the CN production rate as the 
comet approaches the Sun and decreasing faster after perihelion. 
Despite the evolution of the C 2 to CN ratio, comet Lovejoy can 
be classified as a typical comet in terms of carbon-chain species 
as defined by A’Hearn et al. (1995). However, around 2.5 au, 
the C 2 -to-CN ratio almost reaches the lower limit for a typi¬ 
cal comet. This shows that the composition of the coma evolves 
during the orbit and that a comet could be classified as typical 
or depleted depending on the heliocentric distance at which it 
is observed. We must then treat this kind of classification with 
caution. 


Table 4. Pre- and post-perihelion fitted power-law slopes for the 
CN, C 3 , C 2 , and A(0)fp to OH ratios heliocentric dependences. 


Species 

r-dependence 

pre-perihelion post-perihelion 

log[Q(CN)/Q(OH)] 

1.00+0.52 

1.06+0.20 

log[Q(C 3 )/Q(OH)] 

1.02+0.69 

0.70+0.30 

log[Q(C 2 )/Q(OH)] 

-0.35+0.44 

-0.21+0.23 

log[Q(A(0)/p)/Q(OHj] 

-0.23+0.36 

-0.48+0.20 


In Fig. 5 we show the evolution of the CN, C 3 , andC 2 pro¬ 
duction rates, and A(0)fp compared to the OH production rate. 
We did not represent the evolution of the NH/OH ratio since si¬ 
multaneous measurements of these two species were scarce. The 
dispersion of the ratios is mainly due to the dispersion of the 
OH production rates caused by the lower signal-to-noise ratio 
in these images. We fit two separated slopes for pre- and post¬ 
perihelion data, which are given in Table 4. All the ratios dis¬ 
play a slight asymmetry about perihelion, post-perihelion values 
are higher than pre-perihelion values owing to the asymmetry 
of the OH production rates about perihelion. The asymmetry is 
stronger for A(0) fp since, as already noticed, the comet is dustier 
post-perihelion. Except for the asymmetry about perihelion, the 
A(0)fp -to-OH and C 2 -t- OH ratios are almost constant with he¬ 
liocentric distance. For these two ratios, the slopes given in Table 
4 are approximately consistent with 0. In contrast, the slopes of 
the CN-to-OH and, to a lesser extent, C 3 -to-OH ratios are differ¬ 
ent from zero. This means that the evolution of C 2 , OH, and the 
dust are correlated, while C 3 and CN are likely not, as suggested 
by the analysis of the heliocentric dependences in section 3.1. 
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1 1.5 2 2.5 3 

r [au] 


Fig. 4. Ratio of C 2 -to-CN production rates as a function of the 
heliocentric distance. Pre-perihelion values are represented with 
filled symbols and post-perihelion values with open symbols. 



log r [au] 


Fig. 5. Ratio of CN, C 3 , C 2 production rates, and the A(0)fp to 
the OH production rate as a function of the heliocentric distance. 
Pre-perihelion values are represented with filled symbols and 
post-perihelion values with open symbols. Full lines represent 
linear fits of the ratios variation with the heliocentric distance 
pre-perihelion. Dotted lines are the same for post-perihelion val¬ 
ues. 

3.3. Morphology 

In this section, we describe the dust and gas coma morphologies 
of comet Lovejoy. We tested several techniques to improve the 
contrast between the coma and the specific features: subtraction 
and division of an azimuthal median profile, rotational filter, and 
Larson-Sekanina filter (see Schleicher & Farnham 2004; Larson 
& Sekanina 1984). All these techniques give approximately the 
same results, but the removal of the azimuthal median profile 
allows fainter large scale structures to be detected, and we thus 
present the images obtained with this technique. 


We first discuss the observations obtained before perihelion 
between September 12 and November 16, 2013. We detected 
features in all the filters, indicating that the nucleus of comet 
Lovejoy has several active regions. The first CN feature was de¬ 
tected around September 21, 2013, while the dust, C 3 , and C 2 , 
features were only detected later, around October 6 , 2013. The 
OH and NH features do not appear in all the images, owing to 
the high airmass of some observations and the lower signal-to- 
noise ratio in these bands. We present in Fig. 6 the results of 
the subtraction of the azimuthal median profile for C 2 , C 3 , CN, 
NH, and OH images, and also for the dust in RC filter for im¬ 
ages taken on November 2 and 3, 2013. The RC image shows an 
enhancement of the coma in the tail direction around PA 280°. 
There is a broad and slightly curved feature centered on PA 100°. 
The C 2 image shows two jets around PA 100° and 280°, at the 
same position as the dust features mentioned above. We also see 
a broad feature around PA 15° and a fainter one around PA 190° 
that are not visible in the dust images. The features around PA 
100 ° and 280° visible in both C 2 and dust images seem to indi¬ 
cate that at least part of the C 2 originates in the same region(s) 
as the dust, or maybe that the C 2 is released from dust grains. 
The CN morphology is different from the C 2 morphology. CN 
displays a roughly hourglass shape with two main fans around 
PA 15° and 175°. Features were also visible in the C 2 images at 
approximately the same position as the CN jets but were much 
weaker. The C 3 morphology is rather similar to the C 2 morphol¬ 
ogy, with one feature around PA 280° and a large fan around PA 
5°. However, the C 2 feature around PA 100° and the weak one 
around PA 190° are not observed in the C 3 image. To summa¬ 
rize, C 2 and, to a lesser extent, C 3 images are the sum of dust 
and CN features. In the NH and OH filters, we do not detect jets 
but only an enhancement of the coma centered on the tail direc¬ 
tion, around 280°. Since the comet was only visible one or two 
hours at the end of the night during this period, we could not de¬ 
tect variations in the feature’s shape and position during the night 
caused by the nucleus rotation. The position and the shape of the 
features described above for November 2 and 3 varied slightly 
between late September and mid-November because of changes 
in the viewing geometry. We also have images of the comet post¬ 
perihelion. Unfortunately, we could only identify clear features 
for CN. The morphology is almost the same as in November 
2013: an hourglass shape with two main fans around PA 15° and 
205°. 

4. Discussion 

The study of the evolution of comet Lovejoy activity with the he¬ 
liocentric distance revealed a significant asymmetry about peri¬ 
helion. This asymmetry is visible from C 2 , C 3 , CN, NH, and OH 
production rates but also from the dust, even if the amplitude 
of the asymmetry varies depending on the species. The A(0)fp 
values are at least a factor two higher after perihelion than be¬ 
fore at the same heliocentric distance. Asymmetries of the activ¬ 
ity about perihelion seem to be a common phenomenon among 
comets (see, for example, A’Hearn et al. 1995, Schleicher et al. 
1998, Knight & Schleicher 2013, Opitom et al. 2015). It remains 
difficult to know if such asymmetries happen for all comets or 
to understand their origin, because large datasets with sufficient 
temporal coverage and with observations on both sides of per¬ 
ihelion are still too scarce. A possible cause often mentioned 
is a thermal inertia effect. One other probable cause is a sea¬ 
sonal effect on a nucleus with several active regions. Indeed, for 
comets on highly eccentric orbits, changes in the nucleus illu¬ 
mination can happen close to perihelion. If the illuminated per- 
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Fig. 6. C/2013 R1 (Lovejoy) RC, C 2 , C 3 CN, NH, and OH images from November 2 and 3, 2013 processed by subtracting an 
azimuthal median profile. All images are oriented with north up and east left. The field of view is 4.3’ x 4.3’. The arrow indicates 
the direction of the Sun. 
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centage of the active region(s) rises, this will also impact the evo¬ 
lution of the comet gas and dust production rates. Some regions 
might also need more heat to be activated, for instance if they are 
buried deeper, or there might be some delay before they start to 
outgas. The study of comet Lovejoy coma morphology revealed 
gas and dust jets coming out of the nucleus. This strongly sug¬ 
gests the existence of one or several active region(s) on the nu¬ 
cleus. This is consistent with a seasonal effect being the cause of 
the observed asymmetry. We also noted that the dust-to-gas ratio 
changed; the comet is dustier after perihelion than before. This 
may bring us important clues about the heterogeneity level of 
the nucleus. More precisely, this could be explained by different 
ice and dust mixing in the active regions. However, we must be 
cautious while drawing conclusions since we do not completely 
understand the mechanism that causes such asymmetries. 

While looking at the evolution of the comet activity with the 
heliocentric distance, we noticed a correlation between some of 
the radicals and the dust. The slopes of the C 2 , NH, OH, and the 
dust heliocentric dependence presented in Table 3 are consistent 
with each other within the error bars, while the CN slope is shal¬ 
lower. The C 3 slope is close to the CN slope. This correlation 
between some radicals and the dust is also observed after peri¬ 
helion. It is even more obvious looking at the evolution of the C 2 , 
C 3 , CN, NH, and dust production rates relative to OH. Except for 
the asymmetry about perihelion, the C 2 /OH and A(0)/p/OH ra¬ 
tios do not significantly vary with heliocentric distance. On the 
other hand, a power-law fit of the CN/OH and the C 3 /OH ratio 
with the heliocentric distance gives a slope different from zero. 
Once again, these trends are identical on both sides of perihelion. 

All this evidence suggests that the radicals we observe in the 
coma of comet Lovejoy may not all come from the same source. 
The issue of gas coming from icy or carbonaceous grains in the 
coma of comets has been debated for a long time. C 2 H 2 , C 2 H 6 , 
and HCN are usually thought to be the primary parents of C 2 
and CN (Helbert et al. 2005, Weiler 2012). However, discrepan¬ 
cies between the parent and daughter production rates observed 
for some comets seem to indicate that we need other parents to 
explain the abundances of some radicals in the coma (Fray et al. 
2005; Dello Russo et al. 2009; Kobayashi et al. 2010; McKay 
et al. 2014). It has been suggested that C 2 and/or CN can come 
-at least partially- from carbonaceous grains or another unknown 
source. In the case of comet Lovejoy, we see a strong correlation 
between the evolution of the dust production and the evolution 
of some radicals such as C 2 , NH, and OH. This suggests that at 
least part of the C 2 , on one hand, and also part of the NH and 
OH, on the other hand, could be produced by organic-rich and 
icy grains, respectively. Since the evolution of the CN produc¬ 
tion is different from the C 2 and does not seem correlated with 
the dust, most CN would be produced directly by the photodis¬ 
sociation of HCN. 

This hypothesis is further reinforced by the study of the coma 
morphology. Indeed, Fig. 6 shows that the dust and C 2 mor¬ 
phologies have some similarities. C 2 is enhanced in the tail di¬ 
rection, and we see a jet in the sunward direction at the same 
position as the dust jet. We also notice features around the same 
position as the CN fans. In contrast, the CN morphology is com¬ 
pletely different from the dust. It has an hourglass shape with 
two jets around PA 15° and 175°. The C 3 morphology is similar 
to the C 2 morphology. Since HCN, C 2 H 2 , and CofL, usually have 
the same spatial distribution (Dello Russo et al. 2009, Villanueva 
et al. 2011), this could indicate that, while most of the CN is the 
product of HCN photodissociation, only part of the C 2 is pro¬ 
duced by the photodissociation of C 2 H 2 (and C 2 H 6 but C 2 H 6 is 
not very efficient to release C 2 , see Weiler 2012). The rest of the 


C 2 is probably coming from organic-rich grains. This would ex¬ 
plain the enhanced C 2 regions around PA 15° and 175,° while 
the most contrasted features are similar to the dust. We men¬ 
tioned earlier that NH and OH parents could be released from 
icy grains. If this is the case, we would expect the NH and OH 
coma to be enhanced in the tail direction since the icy grains are 
pushed away from the Sun by the radiation pressure. This is ex¬ 
actly what we see in Fig. 6. Both NH and OH are enhanced in the 
anti-solar direction, strongly reinforcing the hypothesis of water 
and ammonia (the most likely parent of NH) being released by 
icy grains. 

To better understand the relationship between parent and 
daughter species, we can compare measurements of the parent 
species production rates made in the IR and the radio domains to 
the radicals production rates we derived here. Only a few results 
have been published so far about comet C/2013 R1 (Lovejoy). 
We summarize the results reported in Paganini et al. (2014b), 
Agundez et al. (2014), and Biver et al. (2014) in Table 5, along 
with our daughter species production rates. From this table, we 
notice that there is a significant dispersion between production 
rates measured almost simultaneously with different instruments 
and from different wavelength ranges. This dispersion is proba¬ 
bly due in part to the different fields of view (less than 10 arc- 
seconds in the millimeter range and 24 arcseconds in the IR) and 
models used but also maybe to the rotational variation. 

We first compare the OH and H 2 O production rates. Our val¬ 
ues agree within the error bars with the ones reported by Biver 
et al. (2014) and Paganini et al. (2014b) for early November. 
In October, however, the OH production rates we derived are 
slightly higher than the water production rates measured by 
Paganini et al. (2014a). 

Since measurements of the HCN production rate of comet 
Lovejoy have also been published by Paganini et al. (2014b) and 
Agundez et al. (2014), the comparison with our CN measure¬ 
ments is of great interest for determining the origin of the CN 
in the coma of comet Lovejoy. In October, the CN and HCN 
production rates are consistent between each other (except for 
October 27 and 29 for which HCN production rates reported by 
Paganini et al. (2014b) are almost a factor two lower than those 
measured a few days before and a few days after, maybe because 
of rotational variation very close to the nucleus). In November, 
our CN production rates agree with those reported by Paganini 
et al. (2014b) and Biver et al. (2014), especially if we take the 
rise in comet activity into account. The HCN production rates 
reported by Agundez et al. (2014) are lower than our CN produc¬ 
tion rates, but they are also lower than those derived by Paganini 
et al. (2014b) and Biver et al. (2014) only a few days earlier. 

Finally, we can also compare the C 2 production rates to the 
production rates of two possible parents, CjHg and C 2 H 2 , re¬ 
ported by Paganini et al. (2014b). For all dates except one, the 
C 2 H 2 production rate upper limit is lower than the C 2 produc¬ 
tion rate. For some dates the C 2 H 2 production rate upper limit is 
almost a factor ten lower than the C 2 production rate. The C 2 H 6 
is of the same order of magnitude as the C 2 production rate. 

The comparison of parent and daughter species production 
rates seems to confirm our interpretation of the origin of daugh¬ 
ter species in the coma of comet Lovejoy. Most of the HCN pro¬ 
duction rates reported are of the same order of magnitude as our 
CN production rates, the ones reported by Paganini et al. (2014b) 
are even in good agreement with our CN production rates. This 
suggests that HCN is probably the main parent of CN in the 
coma of comet Lovejoy. For C 2 , there are large discrepancies be¬ 
tween the most likely parent abundance and the C 2 abundance. 
Indeed most of the upper limits of the C 2 H 2 production rates are 
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Table 5. Comparison of mother and daughter species production rates. 


Date 

Q(OH) a 

Q(H 2 0) B 

Q(CN) a 

Q(HCN) b 

Q(C 2 ) a 

Q(C 2 H 2 ) b 

Q(C 2 H 6 ) b 

2013-Oct-21 

2013-Oct-22 

2380±250 

1850±190 

6.38±0.19 

5.9±1.2 

7.63±0.21 

<9.2 

13.80±1.70 

2013-Oct-24 

2720±290 

1530±250 

6.96±0.21 


9.10±0.34 



2013-Oct-25 

2013-Oct-26 

2860±330 

2320±330 

7.20±0.31 

6 .8±0.6 

9.62±0.39 

<1.5 

21.70±3.50 

2013-Oct-27 


1560±250 


3.6±0.1 


<2.2 

17.40±3.70 

2013-Oct-29 

2013-Nov-2 

2013-Nov-3 

3010±350 

1450±410 

7.97±0.27 

4.7±0.6 

12.20±0.30 

<1.7 

10.80±2.50 

2013-Nov-7 


3680±260 


8.7±0.7 


<1.6 

21.80±2.50 

2013-Nov-8/12 


d 5000 


ct 7.4±0.1 




2013-Nov-13 

3024±790 


10.10±0.70 

c 6.0 

16.10±0.50 



2013-Nov-16 



12.90±1.20 

c 6.1 

20.60±0.90 




All production rates are expressed in units of 10 25 mol/s. 
a This paper 

b Paganini et al. (2014b), from observations with the Near InfraRed SPECtrometer at the Keck observatory 
c Agundez et al. (2014), from Nancay OH observations 
d Biver et al. (2014), from observations with the IRAM 30m telescope 


much lower than the CA production rates. This indicates that we 
need another parent to explain the C 2 abundance in the coma or 
that C 2 comes in part from organic-rich grains. Another possi¬ 
ble parent of C 2 is C 2 H6. Paganini et al. (2014b) detected C 2 H 6 
in the coma of comet Lovejoy, and its abundance is of the same 
order of magnitude as the C 2 . However, it has been reported that 
the release of C 2 from C 2 H6 is not really efficient (Weiler 2012) 
and makes it a poor candidate for C 2 parentage. 

Nevertheless, we must remain cautious when comparing par¬ 
ent molecules production rates measured from the infrared to 
daughter molecules production rates from the optical. The mod¬ 
els used and the parameters in these models are different, making 
any direct comparison difficult. As an example, in this paper we 
used an outflow velocity of 1 km/s, as in A’Hearn et al. (1995). 
However, in most infrared and radio studies, including those 
from Paganini et al. (2014b) and Biver et al. (2014) mentioned 
earlier, an outflow velocity of v = 0.8/' 05 km/s is assumed. To 
assess the importance of the outflow velocity used in the model 
on the derived daughter molecules production rates and then on 
the comparison between daughter and parent species production, 
we computed the gas production rates presented in Table 5 again 
using the expansion velocity law: v = 0.8r 05 km/s. These re¬ 
sults are shown in Table 6. A comparison of Tables 5 and 6 shows 
that the impact of the outflow velocity is significant. The produc¬ 
tion rates derived from optical data are lower if we use the same 
velocity law as for the infrared data (lower by a factor 1.8 in 
October and by a factor 1.5 in mid-November). However, the 
agreement between the observed profile and the Haser model is 
much better for constant expansion velocities of 1 km/s. This 
is not surprising since the Haser scalelengths we use (A’Hearn 
et al. 1995) have been fitted using an expansion velocity of 1 
km/s. 

Examination of Table 6 shows that given the variability of 
water production rates measurements, there are no significant 
differences between OH production rates and those of H 2 O for 
most of the dates. Some of the CN production rates measured 
are lower than those of HCN reported in the literature, but they 
are consistent with the lower HCN production rates measured 
in late October and mid-November. In early November, C 2 pro¬ 
duction rates are still almost a factor five higher than the C 2 H 2 
production rates upper limits. 


Another question is the dependence of both scalelengths and 
fluorescence efficiencies on the heliocentric distance. This of 
course influences the daughter production rates we derive and 
then the comparison between daughter and parent molecules 
production. In our case, we only compare daughter and parent 
molecules for observations made between 1.4 and 1 au from the 
Sun so that variations from the r 2 scaling used for both scale- 
lengths and g factors should not have a strong influence on our 
results. For example, if we use a r 3 scaling law for the scale- 
lengths, the maximum variation in the production rates over this 
range heliocentric distances is 20%. Finally, the C 2 fluorescence 
efficiencies we used only consider the C 2 in triplet state since 
the ratio of C 2 in triplet and singlet state and its heliocentric de¬ 
pendence are not understood well. We thus underestimate the C 2 
total production rates, which reinforces our conclusion about C 2 
parentage. 

5. Summary 

We performed a photometry and imaging monitoring of the Oort 
Cloud comet C/2013 R1 (Lovejoy) during several months with 
the robotic TRAPPIST telescope at La Silla observatory. We ob¬ 
served the comet on both sides of perihelion with good temporal 
coverage. We derived production rates and Afp values and stud¬ 
ied the evolution of its activity. We also performed a morpholog¬ 
ical analysis of the coma. 

The morphological analysis revealed several features in the 
coma, both for the dust and the gas species. This indicates that 
the nucleus has one or several active region(s). The comparison 
of the comet activity pre- and -post perihelion shows an asym¬ 
metry of the gas and dust production rates about perihelion, the 
rate of brightening is steeper than the rate of fading. We esti¬ 
mated that it is likely to be seasonal effect caused by changing 
system geometry and insolation of the active region(s). The dust 
production rate is at least two times higher after perihelion than 
before at the same heliocentric distance, which suggests a certain 
degree of heterogeneity of the nucleus. 

The evolution of the comet activity and production rates ra¬ 
tios with the heliocentric distance, the comparison of parent and 
daughter species production rates, as well as the coma mor¬ 
phology all indicate that the C 2 , NH, and OH productions are 
strongly correlated to the dust production, while the CN pro- 
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Table 6. Comparison of mother and daughter species production rates using a gas outflow velocity 
t' = 0.8r -0 - 5 km/s to compute daughter species production rates. 


Date 

Q(OH) a 

Q(H 2 0) b 

Q(CN) a 

Q(HCN) b 

Q(C 2 ) a 

Q(C 2 H 2 ) b 

Q(C 2 H 6 ) b 

2013-Oct-21 

2013-Oct-22 

1280±130 

1850±190 

3.53±0.11 

5.9±1.2 

4.20±0.12 

<9.2 

13.80±1.70 

2013-Oct-24 

1500±160 

1530±250 

3.95±0.12 


5.15±0.19 



2013-Oct-25 

2013-Oct-26 

1610± 190 

2320±330 

4.16±0.18 

6 .8±0.6 

5.52±0.22 

<1.5 

21.70±3.50 

2013-Oct-27 


1560±250 


3.6±0.1 


<2.2 

17.40±3.70 

2013-Oct-29 

2013-Nov-2 

2013-Nov-3 

1810±210 

1450±410 

4.90±0.17 

4.7±0.6 

7.46±0.16 

<1.7 

10.80±2.50 

2013-Nov-7 


3680±260 


8.7±0.7 


<1.6 

21.80±2.50 

2013-Nov-8/12 


cl 5000 


■^AtO.l 




2013-Nov-13 

2090±510 


6.88±0.44 

c 6.0 

10.90±0.40 



2013-Nov-16 



8.97±0.84 

c 6.1 

14.30±0.60 




All production rates are expressed in units of 10 25 mol/s. 
a This paper 

b Paganini et al. (2014b), from observations with the Near InfraRed SPECtrometer at the Keck observatory 
c Agundez et al. (2014), from Nancay OH observations 
d Biver et al. (2014), from observations with the IRAM 30m telescope 


duction is not. The case of C 3 is more complicated and quite 
intermediate between these two behaviors. We then suggested 
that an important part of the C 2 could be released from organic- 
rich grains and that water and ammonia (the parents of OH and 
NH) could be released from icy grains. CN and some C 2 would 
be produced by the photodissociation of HCN and C 2 H 2 , which 
usually display the same spatial distribution. C 3 would also be 
produced partly from the organic-rich grains and by the pho¬ 
todissociation of a parent molecule. 
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